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The  transient  response  of  a system  of  two  initially  concentric  circular 
cylindrical  elastic  shells  coupled  by  an  ideal  fluid  and  l^ilnged  by  an  incident 
plane  pressure  pulse  is  stuped.  The  classical  techniques  of  separation  of 
variable  and  Laplace  transforms  are  employed  for  simultaneously  solving  the 
wave  equations  governing  the  fluid  motions  and  the  shell  equations  of  motion. 

The  transformed  solutions  are  arranged  in  such  a manner  that  their  inverse 
transforms  can  be  accurately  calculated  by  solving  a set  of  Volterra  Integral 

DO  t jan  *?i  1473  EOlTlON  OF  » MOV  At  It  OBSOLETE  \ _ 

S N 0102-014- AAOI  _ 

SECURITY  CLASSIFICATION  of  THIS  PAGE  (Whan  Data  Bnt  at  ad) 


$,£%  ISO 

i 1 0 


-WWlT'f  CLASSIFICATION  OF  THIS  PAOEWw  Dsfs  Bntmrmd) 


Block  20 


equations  In  the  time  domain.  A sample  calculation  of  shell  response  was 
performed  and  results  are  compared  to  the  case  In  which  the  outer  shell  Is 
absent.  It  Is  found  that  the  primary  effects  of  a thin  outer  shell  could  be 
estimated  by  simple  asymptotic  formulae. 
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TRANSIENT  RESPONSE  OF  TWO  FLUID-COUPLED  CYLINDRICAL  ELASTIC 
SHELLS  TO  AN  INCIDENT  PRESSURE  PULSE 

INTRODUCTION 

In  order  to  gain  physical  Insight  In  the  response  as  well  as  to 
provide  a data  base  for  the  development  of  general  numerical  methods  for 
predicting  the  underwater  explosion  response  of  fluid  coupled  shell 
systems  such  as  the  double  hull  section  of  a submarine,  the  transient 
response  of  systems  of  fluid  coupled  elastic  shells  of  simple  con- 
figurations to  the  excitations  of  incident  pressure  pulse  are  being 
analytically  investigated.  The  case  of  two  concentric  spherical 
elastic  shells  lends  itself  to  the  classical  treatment  of  the 
separatlon-of-variable  and  integral  transform  techniques  and  accurate 
solutions  have  been  obtained  by  a successive  integration  scheme  [1]. 

It  was  found  that  a thin  outer  shell  tends  to  be  transparent  to  the 
incident  pulse. 

In  the  present  report,  a solution  method  and  results  for  the  case 
of  two  concentric  cylindrical  elastic  shells  of  infinite  length  are 
presented.  This  is  also  one  of  the  cases  in  which  it  is  possible  to 
apply  the  separation-of-variable  and  Laplace  transform  techniques  to 
simultaneously  solve  the  wave  equations  governing  the  wave  motions  in 
the  fluids  and  the  eauations  of  motion  of  the  elastic  shells.  The 
corresponding  problem  for  a single  cylindrical  shell  has  been  widely 
studied  and  the  results  have  been  reviewed  and  summarized  in  many 
articles,  e.g.,  references  [2,3].  Here,  it  was  more  convenient  to 
obtain  the  inverse  Laplace  transform  for  the  transient  solution 
indirectly  using  a integral  equation  method  [4]  or  a dif ferential- 
Manuacript  aubmitted  June  30,  1978. 
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integral  equation  method  [5].  It  is  also  possible  to  formulate  the 
present  double-cylindrical-shell  problem  such  that  the  same  integral 
equation  method  can  be  used  to  calculate  the  transient  shell  response. 

Results  are  obtained  for  a system  of  two  steel  shells  submerged 
in  water.  Similar  to  the  spherical  case  [1],  it  is  also  found  that  a 
thin  outer  shell  does  not  change  significantly  the  basic  response 
feature  of  the  inner  shell  and  that  its  primary  effects  are  to  reduce 
the  axi-synmetric  (breathing)  and  translational  motions  of  the  inner 
shell.  These  lower  mode  effects  can  be  estimated  by  simple  formulae. 
DESCRIPTION  OF  THE  PROBLEM 

Figure  1 sketches  the  submerged  fluid-coupled  cylindrical  shell 
system  and  the  incident  plane  pressure  wave.  The  fluid  surrounding  the 
outer  shell  and  that  between  the  two  shells  are  considered  to  be  ideal 
compressible  fluids  in  linear  wave  motions  and  can  be  characterized  by 
their  unperturbed  mass  densities  and  sound  speeds  (p  , c ) and  (p , c) 
respectively.  The  shells  are  initially  concentric.  In  this  study, 
the  strength  of  the  incident  wave  is  sufficiently  weak  such  that  the 
shell  deflections  are  elastic  and  small  and  deviation  from  the  con- 
centricity remains  negligible  for  the  time  duration  of  interest.  The 
middle  surface  radii,  thickness,  mass  densities,  Young's  moduli  and 
polsson's  ratios  of  the  outer  and  inner  shells  are  (a  , h , pg,  E , v ) 

and  (a,h,p  ,E,v),  respectively.  The  z-coordinate  of  the  cylindrical 

s 

coordinate  system  (r,  6,  z)  coincides  with  the  axis  of  the  shells  in  its 
unperturbed  position.  The  incident  plane  wave  front  is  also  parallel  to 
the  z-axis  and  therefore  all  response  entities  do  not  vary  with  respect 
to  z. 
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The  deflections  of  the  inner  shell  in  the  r-  and  0-  direction, 


e 

normalized  with  respect  to  the  outer  shell  radius  a , are  denoted  by  w 

6 6 

and  v respectively  and  those  of  the  outer  shell  by  w and  v respectively. 
The  deflections  in  the  z-direction  do  not  enter  into  the  present  two- 
dimensional  system  of  equations  of  motion.  The  total  pressure  field 
exterior  to  the  outer  shell  is  denoted  by  pe(r,  0,  t)  and  that  between 
the  shells  by  p (r,  0,  t)  where  t designates  time.  The  following 
dimensionless  parameters  are  also  used  in  the  mathematical  formulation: 


and 


respectively,  where  V2  is  the  Laplacian  operator.  The  boundary 
conditions  of  the  problem  are  that  P6  satisfies  the  radiation  condition 
at  far  field  and  that 


3P  _ 3Pf_ 
9R  3R 


at  R = 1 


3P 

■Ml  “ P. 


at  R = C . 


All  quantities  except  the  incident  pressi  field  have  quiescent 
initial  conditions. 


A Laplace  transform  pair  is  defined  as 


ao 

w (e,s)  = / w (0 ,T)  e sT  d 

A)  T (6) 

y +i°° 

w (0,T)  =/  w (0,s)  e ds 

y-i00 

where  y lies  to  the  right  of  all  singularities  of  w in  the  complex  s- 
plane  and  i = (-1)^ 

Due  to  the  circular  geomtry  of  the  problem,  the  solutions  can  be 
expanded  in  Fourier  series  as  the  following: 


P(R,0,T)  = £ P (R.T)cos  n 0 
n=0  n 


Pe(R,0,T)  = I P6  (R,T)cos  n 0 
n 

n=® 


w (0,T)  = £ w (T)cos  n 0 
n-0  n 


v (6,T)  = £ v (T)sin  n 0 
n-1 

CO 

we  (Q.T)  - S we  (T)coa  n 9 
n“0  n 
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V (0,T)  - E V ' (T)sin  n 0 
i n 
n-1 


This  type  of  expansion  is  not  suitable  for  the  pressures  and  shell 
accelerations  at  early  time  [3,4].  Otherwise,  it  greatly  facilitates 
the  calculation  of  the  responses  of  the  shells. 

For  the  solution  scheme  to  be  used  here , any  linear  elastic  theory 
of  the  cylindrical  shell  is  applicable.  In  order  to  juxtapose  solutions 
to  those  for  a single  shell  previously  obtained  in  Ref.  [4],  the  version 
of  shell  equations  of  motion  used  therein  is  also  used  here.  In  the 


Laplace  transform  domain,  they  are 
-Ms (s2s  + a)  P (S ,s) 


S4  s4  + + S2s2  + p 
n n 


(8) 


-Me(s2  + a*)  l?J  (l,s)  - Pn(l,s] 

S4  + +e  S2  + pC 

n n 


and 


n - 0,1,2,... 


(S2s2+a  ) v - -n  C2  (l+n2I)  w 
n n n 


(9) 


(s2  + a6)  v6  ■ -n  Cz  (l+n2Ic)  wc 
n n e n 


’2  /1-W.2-re^ 


n - 1,  2,  3,  ... 

The  total  pressure  field  exterior  to  the  outer  shell  consists  of 
the  pressure  due  to  the  incident  wave  and  those  due  to  scattering  and 
radiation  by  the  outer  shell.  An  arbitrary  incident  plane  pressure 
wave  impinging  the  outer  shell  at  (0*0,  R“l)  at  T“0  can  be  expressed 
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by  the  following  aeries  [4]. 


P*  (R,0,s)  = f(s)e  8 Z_  e I (sR)cos  n 0 

n=0  n n 


(10) 


where 


f(8>  ia  the  Laplace  transform  of  the  temporal  characteristics  of  the 

incident  wave,  I (sR)  is  the  modified  Bessel  function  of  th <-  first  kind 
n 

and 


c = 1 for  n = 0 
n 


2 for  n > 0 . 


(11) 


SOLUTIONS  IN  THE  LAPLACE  TRANSFORM  DOMAIN 

It  can  be  readily  shown  that  the  solutions  to  the  system  of  Eqs. 
(2)  through  (11)  are: 


P (R,s)  = f (s)  e e 
n n 


VR> 


I (s)  K (sR)  1 s^ 
n n n 

Ms)  -I  K (s)  n 


(sR) 


(12) 


and 


sw  K (c  £s)  - sw  K (c  s;  / \ T , nv 

F (R,s)  ■ ° ° T - — Q — E — / ££_  \ WR) 

n ’ » / s I e el 

An  <8>  \P  c / 

~ — e j v 

sw  I (c  s)  - sw  I (c  5s)  / \ 

+ — n ■" I n n.  - -r  — [ ) K (c  8 R) 

lee/  n ' r 
A (s)  \ P c / 


w®  - M®(s2-kx®)  [(^4s4++  S2s2+u  ) A (s)(p  c )/(pc) 

n n n n n 

f(s)e-8e 

- M<;s(c282+an)Bn(s)]  A (s— 

n 


f(s)e”Se 

wn  - MMe(c2824an)(82+an) 

r n 


(13) 


(14) 


(15) 


--  ■ r -• 


where  K (s)  is  the  modified  Bessel  function  of  the  second  kind,  the 
n 

prime  denotes  differentiation  of  the  Bessel  functions  with  respect  to 
their  arguments, 


A (s)=s (£4s4+t  S2 s2+y  ) [K  (s)  (s‘4++es2+pe)-sK  (s)^ (s2+«e)  ] A (s)  (p6ce)  / (pc) 
n nnn  nnn  nn 

-s2K  (s)  [Me (s2+a*')  (C4s4-H  C2s2+y  )G  (s)+M5(;2s2-hx  ) (s4++®s2-hj®) B (s)  ] 
n n nnn  nnnn 

+MMecs3(;2s2-fa  ) (s2+a®)  [K  (s)B  (s)+K  (s)L  (s)  (pc)/(pCce)  ] (16) 

n nnn  n n 

A (s)-I  (c  Cs)K  (c  s)-I  (c  s)K  (c  ts)  (17) 

n nr  nr  nrnr 

B (s)-I  (c  CS)K  (c  s)-I  (c  s)K  (c  cs)  (18) 

n nr  nr  nrnr 

G (s)-I  (c  s)K  (c  Cs)-I  (c  Cs)K  (c  8)  (19) 

n nrnr  nr  nr 

and 

L (s)«I  (c  s)K  (c  Cs) -I  (c  Cs)K  (c  s)  . (20) 

n nrnr  nr  nr 


v and  v can  be  found  from  Eq.  (9).  Strains  and  stresses  of  the 
shells  can  be  computed  using  appropriate  strain-deflection  and  stress- 
strain  relationships. 

With  use  of  the  Tauber's  theorem  of  Laplace  transform  [16],  some 
of  the  asymptotic  behaviors  of  the  shell  responses  at  late  time  can 
readily  be  revealed  from  Eqs.  (14)  and  (15).  Specifically,  for  the 
case  where  the  incident  wave  P*  is  a unit  step  wave,  i.e. , f(s)  ■ l/s. 


-M6 [c2C2 (l-c2)+2pr£2M] 
c2C2C2(l-C2)+2prMeC2+2pr;2MC2 


(21) 
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?■ 

% 


1.^ 


VT) 


-2pr5MMe 

c2c2c2(l-^2)+2prMeC2+2pr?2MC^ 


-4Me[(l-C2)-K).5p  M(l+42)] 

w (T)  = 

1 2 (1— C2 ) (2+Me)+2p  (1+42) (M+Me)4p  MMe[l+42+c  (1-42)] 

X -»«»  r r r 


wx(T)  = 


-4p  MM 
r 


T •*>  2(l-42)(2+Me)+2pr(l+42)(M«ie)+PrMMe[l+42+crCl-C2)] 


where  the  dot  denotes  differentiation  with  respect  to  T. 

Eqs.  (21)  and  (22)  give  the  late  time  shell  deflections  long  after 
the  incident  wave  has  engulfed  the  outer  shell  and  they  can  also  be 
obtained  by  static  analyses.  Eqs.  (23)  and  (24)  are  expressions  for  the 
late  time  translational  velocities  of  the  shells  in  the  direction  of 
propagation  of  the  incident  wave.  It  can  be  shown  that  the  only 
condition  under  which  w^(<*>)  = w®  (“)  is  PrM=2,  i.e.  , when  the  inner 
shell  is  neutrally  buoyant  in  the  interior  fluid.  Otherwise,  the 
ratio  between  them  depends  only  on  the  buoyancy  of  the  inner  shell  and 
the  radius  ratio  4.  In  Fig.  2,  it  can  be  seen  that  a positively 
buoyant  inner  shell  (p^M  > 2)  translates  faster  and  a negatively 
buoyant  one  (p^M  < 2)  translates  slower  than  the  outer  shell  at  late 


time . 


Formulae  (21)  through  (24)  can  be  used  to  facilitate  parametric 


studies  of  the  effects  of  various  shell  and  fluid  properties  on  the 
axisymmetric  and  translational  responses  of  the  shells.  They  also 


provide  asymptotic  checks  for  the  numberical  calculations. 
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THE  INVERSION  OF  THE  TRANSFORMED  SOLUTIONS 


The  transformed  solutions  in  Eqs.  (12)  through  (15)  are  highly 
complicated  functions  of  s.  To  find  their  inverse  transforms,  it  is 
more  convenient  to  rearrange  them  such  that  the  solutions  in  the  time 
domain  can  be  obtained  by  numerically  solving  a Volterra  type  of 
integral  equation  the  kernel  of  which  can  be  accurately  calculated. 
This  procedure  has  been  demonstrated  to  be  quite  effective  in  previous 
studies  of  the  single  cylindrical  shell  problem  [4,5]. 

To  proceed,  Eq.  (15)  is  rearranged  to  assume  the  following 


form: 


w (s)  - w (s) 
n n 


Me8(s2+C6)  K (s) 


n n 


(s4++es2+pe)K  (s) 
n pn'  n 


1 -l 


/pc  \M^(^2s\)  V8) 

\pece/(c4s4++  c2s2+v  )An^ 


+ /_£C_\  ^^2«2-*n>  V*>  /pc\  “ 

ypec6/  (C4s4+fnC2s2+vin)  An^  \pEce/  ( 


Mes(s2+a  ) G (s) 


a,  ,e  A (®) 

sH++  s^+p  ) n 
n n 


2 MCs(G2s2+a  ) Mes(s2+a®)  L (s) 

/ pc  \ n n'  n 

W) 


U4s4++  ;2s2+p  ) (s4++es24ve)  An(8) 

n n n n 


p M(c2s2-h*  ) Me(s2+a®) 

r n n 


(C4s4+t  C2s2+u  ) (s4++®s2+v®)  c2A  (s) 
n n n n r n 


f(s)e  8e 


8 K (s) 
n 


(25) 


The  Nielsen's  function, 


Wn(l,T) 


j_  rH' 

2,1  i , 

y-ioo 


Kn(8) 


+ 1 


lVs> 


,T. 

e ds  , 


(26) 


has  been  accurately  computed  and  tabulated  by  Nielsen  [7],  The 
resultant  of  the  incident  and  scattering  transient  pressure  acting  on 
a rigid  and  motionless  circular  cylinder  Impinged  upon  by  an  incident 


4 


t 


r 


i 

i 

' 

| 


plane  wave, 

,rig 


P“°  (T) 
n 


-J -fr 

2tt1  / 


_ s 

Y+i«>  f(s)e  e 


sK  (8) 
n 


n sT, 
— e as 


(27) 


has  also  been  accurately  computed  with  the  aid  of  the  Nielsen's 
function  [4]. 

In  the  following  inverse  Laplace  transform  integrals, 
/Y+1“  Me8(s2-hxe)  eST 


rn  <T> 


— / 

2ui  / 

J v-ia 


SI*++e82+p 

n n 


ds  , 


(28) 


rS  <« 


1°° 

-A-  f / _pc\ 

21,1 4-i.  W) 


MC8(c2s2+an)  Bn(s)  ^ 

(C4s-+f  C2s2+u  )V8) 

n n' 


(29) 


rn  (T) 


2ni  / 

-/Y-i°° 


M s(s2+a  ) G (s) 


n'  n sT , 

1 — rr  e ds  , 

eel  , 4 , . e 2 . e.  A (s) 
pc  J (84++  8z4-p  ) n 


(30) 


r (T)  + 
n 


f Y+i°°  ; 

— / 

2nl  y-ia>  \Pece/ 


2 McsU2s2+am)  Mes(s2+a®)  Ln(s)  gT 

e ds , 


(C4s4++  52s2+p  ) (s4++es2+ue)^n^S^ 


and 


Fn  <T> 


/Y+i» 

-1. 


p M(t2s2+ct  ) Me(s2+aC)  sT 

r n n e 


(32) 


°°  (5484++  c2s2+p  ) (s4++es2-Hie)  c2A  (s) 

vs  n n n n ' rn 


ds , 


the  Integrands  are  single-valued  and  analytic  in  the  complex  s-plane 

since  it  can  be  shown  that  A (s) , B (s) , G (s)  and  L (s)  defined  in 

n n n n 

Eqs.  (17)  through  (20)  are  single-valued  and  analytic.  Their  poles  are 
determined  by  the  zeros  of  An(s)  and  the  4th  order  polynomials  from  the 
shell  equations  of  motion.  By  the  relationships  between  the  Bessel 

10 


(31) 


» 1 


f 


functions  J and  Y and  the  modified  Bessel  functions  I and  K , 
n n n n’ 


- A (±is)  - J (cts)Y(cs)  - J (cs)Y(cCs) 
itn  nr  nr  nrnr 


(33) 


where  the  cross  product  frequently  occurs  in  applied  problems 
involving  an  annular  region  and  therefore  its  zeros  have  been  studied 
and  calculated  by  many  authors.  Ref-.  [8]  is  one  of  the  many 
publications  presenting  tabulated  tables.  The  zeros  of  this  cross 
product  are  infinite  in  number  and  are  all  real  end  simple,  hence 
An(s)  has  infinite  number  of  conjugated  pairs  of  imaginary  zeros. 

The  integrals  in  Eqs.  (28)  through  (32)  can  now  be  accurately 
evaluated  by  the  customary  method  of  residues  and  they  are  sums  of 
sine  of  cosine  functions  of  T.  Thereupon,  the  inverse  Laplace  trans- 


form of  Eq.  (25)  is  a Volterra  integral  equation  of  the  second  kind, 

T 

(T)  (34) 

u a if 

0 


wn (T)~/  Sn(T-t)wn(t)dx  - Fn< 

A 


where 


S (T)  - A (T)  - Y (T)  + Tb(T)  + TC  (T)  - Td  (T) , 


and 


Fn(T)  -jf  ’ 

VT>  •/’ 

Jo 

t 

VT>  ■ / 

J f\ 


T®  (T-t)  Prl8  (T)dT, 
n n 


r®  (T-t)  w (i,t)dx  - r®  (T), 
n n n 


(35) 

(36) 

(37) 


A (T-x)  r“  (x)dx  . 
n n 
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(38) 


_ 


F (T)  and  S (T)  are  continous  and  well  behaved  functions  of  T.  The 
n n 

numercial  methods  for  solving  Eq.  (34)  have  been  extensively  studied 

0 

and  well  developed  [9],  The  solution  for  w (T)  can  be  likewise 

n 

calculated  and  it  is  also  evident  form  Eq.  (9)  that  v (T)  and  v6(T) 

n n 

can  be  obtained  by  simple  convolutions. 

RESULTS  AND  DISCUSSIONS 

The  integrals  in  Eqs.  (36)  through  (38)  are  evaluated  numerically 

by  the  trapexoidal  rule.  Eq.  (34)  is  solved  using  a simple  linear 

multistep  method  in  which  the  integral  is  also  computed  by  the 

trapezoidal  rule.  Results  are  obtained  for  a case  in  which  both  shells 

are  made  of  steel  and  both  fluids  are  water.  The  material  properties 

and  dimensions  involved  are  such  that 

c - p = 1,  c2  =*  c2  = 12.58121 
r Kr  ’ e 

M6  = 22.09450  M = 4.41890  (39 

he/ae  - 0.00581  h/a  - 0.02905 
C “ 0.8  n =*  0.2 

The  incidence  is  a step  wave  with  f(s)  = 1/s.  For  all  integrations,  a 
time  step  AT  - 0.0125  is  used. 

In  all  subsequent  figures,  the  present  solutions  are  plotted  in 
dotted  lines  and  juxtaposed  to  those  previously  obtained  in  Ref.  [4] 
and  plotted  in  solid  lines  for  the  case  wherein  the  outer  shell  is 
absent. 

Figure  3 shows  the  results  for  the  wq's,  w^'s  and  W2's.  It  can 
be  seen  that  the  numerical  results  for  the  w^'s  and  w^s  correctly 
approach  their  asymptotic  values  as  calculated  by  Eqs  (22)  and  (24) 
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respectively.  At  early  time,  the  presence  of  the  outer  shell  seems 
to  have  little  effect. 

Physically,  the  motion  of  the  outer  shell  starts  at  T * 0 while 

that  of  the  inner  shell  at  T - n ” 0.2.  This  is  also  evident  from  the 

mathematical  properties  of  Eqs.  (17),  (32)  and  (36).  Their  large  time 

translational  velocities  from  Eqs.  (23)  and  (24)  are  respectively 
• • 

w®(»)  - -1.17516  and  w^(<»)  * -1.30361.  If  the  magnitude  of  the 
Incident  pressure  is  such  that  p^[pG(ce)2]  * 1 x 10  \ it  can  be 
estimated  that  at  T-10  the  relative  displacement  between  the  axes  of 
the  two  shells  is  still  one  order  of  magnitude  smaller  than  the 
thickness  of  the  inner  shell.  Therefore,  for  such  an  instance,  the 
deviation  from  concentricity  of  the  outer  and  inner  shells  would  be 
undoubtedly  negligible. 

The  present  differs  slightly  from  that  of  the  single  shell 

case  and  the  presence  of  the  outer  shell  seems  to  have  reduced  the 

rate  of  decreasing  of  at  large  time  from  its  peak  value. 

The  w ' s for  higher  n's  are  plotted  in  Figures  4 and  5.  For 
n 

0 < T 4z  2,  the  present  w^  through  Wy  are  almost  the  same  as  those  for 

the  case  wherein  the  outer  shell  is  absent.  Similar  to  the  single- 

shell  case,  these  wn's  have  markedly  different  oscillatory  characteristics 
* 

for  T > 2 and  their  oscillation  amplitudes  are  decreased  and  the 
periods  increased  by  the  presence  of  the  outer  shell.  Since  the 
magnitudes  of  these  wn's  are  much  smaller  than  those  of  the  lower  n 
terms  to  begin  with,  the  differences  caused  by  the  outer  shell  would 
be  quite  inslgnlcant  for  the  calculation  of  w(6,T).  It  can  also  be 
observed  from  Figures  3,  4 and  5 that  w converges  quite  fast  and  it 


would  suffice  to  use  only  eight  terms  In  the  series  of  Eq.  (7)  to 
calculate  the  transient  displacements,  velocities  and  strains  of  the 
shells  for  all  values  of  6 and  T.  To  calculate  the  early  time 
acceleration  and  pressure,  It  Is  a customary  practice  to  transform 
Eq.  (7)  by  Poisson's  summation  formula  Into  an  Integral  which  Is  then 
asymptotically  evaluated  for  large  8 and  therefore  small  T by  the  saddle 
point  method.  This  has  been  elegantly  carried  out  for  the  single-shell 
case  In  Ref.  [10].  The  application  of  this  procedure  to  the  double- 
shell case  is  rather  complicated  and  is  outside  the  scope  of  the  present 
report. 

Figure  6 shows  the  time  histories  of  the  radial  velocities  at 

various  locations  of  the  inner  shell.  The  presence  of  the  outer  shell 

causes  little  change  of  the  early  time  values  and  the  profiles  of 
• • 

w(6,T).  Other  than  0 ■ tt/2,  the  values  of  w(0,T)  at  later  time  are  re- 
duced primarily  due  to  the  reduction  of  the  w^  term  previously  shown  In 
Figure  3. 

The  hoop  stress  resultant  N(0,T) (force  per  unit  length)at  the 
middle  surface  of  the  Inner  shell  Is  calculated  by 

V6-«  - tV  <* + i?>  <40> 

l-vz 

and  Its  time  histories  at  0 ■ 0,  ir/2  and  it  are  plotted  In  Figure  7. 
Again,  the  shielding  of  the  outer  shell  has  little  effect  on  the 
stresses  at  very  early  time  and  does  not  change  the  profiles  of  the 
stress  time  histories.  At  later  time,  the  stresses  are  reduced  due  to 
the  reduction  of  the  w^  term.  The  present  results  also  clearly 
demostrates  the  generation  of  circumferential  stress  waves  at  the  Inner 
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shell  as  can  be  observed  from  the  repeated  appearance  of  peaks  in  the 
stress  time  histories  at  appropriate  time  intervals.  These  stress 
waves  have  the  same  features  as  those  of  the  single-shell  case  dis- 
cussed in  Reference  [3,4].  They  encircle  the  cylinder  with  the 
dilatational  wave  speed  in  both  directions.  Since  the  attenuation  due 
to  radiation  is  small,  many  circum-navlgatlon  may  be  observed. 

Similar  to  the  double-spherical-shell  case  [1],  it  could  also 
be  summarized  here  that  the  primary  effects  of  a thin  outer  shell  are 
to  reduce  the  breathing  and  translational  motions  of  the  inner  shell 
and  these  could  be  simply  calculated  from  Eqs.  (22)  and  (24). 
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Fig.  2 - Ratios  of  the  late-time  translational  velocities 
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Fig.  6 - Time  histories  of  redial  velocities  of  the  inner  shell 
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Fig.  7 - Time  histories  of  the  hoop  stresses  st  the  middle 
surface  of  the  Inner  shell 
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U S ARMY  WATERWAYS  EXPERIMENT  STN 
P.O.  BOX  631 
VICKSBURG  MS  39180 
ATTN  J STRANGE 
ATTN  W FLATHAU 
ATT.M  TECH  LIB  (UNCL  ONLY) 


AF  INSTITUTE  OF  TECHNOLOGY.  AU 
WRIGHT  PATTERSON  AFB,  OH  464 33 

ATTN  LIB  AFIT  BLDG  640  AREA  BIUNCL  ONLY) 


NKF  ENGINEERING  ASSOCIATES.  INC 
8720  GEORGIA  AVE 
SILVER  SPRING.  MD.  20910 

ATTN  DR.  R.  O.  BELSHEIM 


M&TCO 
2130  ARCH  ST. 
PHILA.  PA.  19103 
ATTN.  B.  McNAIGHT- 
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